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We present the results of the search for rare particles (magnetic monopoles, nuclearites, WIMPs and LIPs) with the MACRO 
detector. For magnetic monopoles (the main goal of the experiment) our limit is ~ 0.4 times the Parker bound for < f3 < 10 _1 . 
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1 Introduction 

MACRO EH at the Gran Sasso Laboratory is a large 
area underground detector devoted to the search for rare 
events in the cosmic radiation. It is optimized to search 
for GUT magnetic monopoles (MM), but can also per- 
form many observations relevant to astrophysics, nuclear, 
particle and cosmic ray physics. The main MACRO 
physics items arc the study of atmospheric neutrinos and 
oscillations, the search for MMs, the study of the high en- 
ergy underground muons, the primary cosmic ray com- 
position, the measurement of the muon residual energy 
spectrum, the search for low energy stellar collapse neu- 
trinos and the high energy neutrino astronomy. 

Here we present the results of the search for rare par- 
ticles: magnetic monopoles and nuclearites, weakly in- 
teracting massive particles and lightly ionizing particles. 
All these searches gave (until now) null results, setting 
significant limits on the fluxes of these rare particles. 

2 The MACRO detector 

The MACRO detector consists of six supermodules (to- 
tal sizes 77 x 12 x 9 m 3 ), each one divided in a lower and 
an upper part. The lower part is made by ten horizontal 
planes of limited streamer tubes, interleaved with seven 
rock absorber layers, and two liquid scintillation counter 
layers on the top and bottom. The lateral walls are closed 
by four "vertical detectors" , formed by a liquid scintilla- 
tor layer sandwiched between two sets of streamer tubes 
(three planes each) . The upper part is made by two "ver- 
tical detectors" on the East and West faces and by a roof 
with one layer of scintillators sandwiched between two 
planes of tubes; this part is left open on the North and 
SouthJaces to house the electronics. A nuclear track de- 
tector! 3 ] is located horizontally in the middle of the lower 
part and vertically on the East and North walls. 

The scintillation counters are equipped with specific 
triggers for rare particles, muons and stellar gravitational 
collapse neutrinos and by 200 MHz WFDs. The streamer 
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tubes are read by 8-channel cards which discriminate the 
signals and send the analog information (time develop- 
ment and total charge) to an ADC/TDC system; the 
discriminated signals form two different chains of TTL 
pulses, which are the inputs for the streamer tube Fast 
and Slow Particle Triggers. 

3 Magnetic Monopoles 

Massive (Mm ~ 10 17 GeV/c 2 ) magnetic monopoL 
arise spontaneously in Grand Unified Theories (GUTs) 
of electroweak and strong interactions. Magnetic 
monopoles of such a large mass cannot be produced with 
accelerators and must be searched in the cosmic radia- 
tion. The MACRO experiment was designed to be sen- 
sitive to. monopoles at a flux level well below the Parker 
Bound I <&m £ 10~ 15 cm~ 2 s -1 sr -1 in the monopole 
velocity range 4 x 10~ 5 < (3 < 1. The use of three subde- 
tectors (liquid scintillators, streamer tubes and nuclear 
track detector) ensure redundancy of information, multi- 
ple cross-checks and independent signatures for possible 
monopole candidates. The results reported here are ob- 
tained using the various subdetectors in a stand-alone 
and in a combined way. All the limits refer to monopoles 
with unit Dirac magnetic charge (g = 137/2 e), catal- 
ysis cross section a < 10 mb (we do not consider the 
monopole induced nucleon decay) and isotropic flux (we 
consider monopoles with enough kinetic energy to tra- 
verse the Earth); the last condition sets a /3-dependent 
mass threshold of ~ 10 17 GeV for f3 ~ 5 x 10~ 5 and lower 
(down to 10 10 GeV) for faster monopoles. 

3.1 Searches with the scintillator subdetector 

The energy loss, arrival time and velocity of a particle 
passing in the scintillator system are measured by using 
the total charge and shape of the photomultiplier pulses. 

In the low velocity region (1.8 x 10~ 4 < (3 < 3 x 10~ 3 ) 
we studied the photomultiplier waveforms, looking for 
the wide, flat and small amplitude signals (or long 
trains of single photoelectrons) expected for a monopole. 
No candidates were found in two independent analy- 
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ses&Q; the flux upper limits (90 % C.L.) are 5.6 and 

4.1 x KT 15 cm" 2 s" 1 sr" 1 (curves "A" and "B" in fig. |l|). 
In the medium velocity range (1.2 x 10~ 3 < (3 < 

1CP 1 ) the monopoles are searched using the data col- 
lected between October 1989 and March 1998_by the 
stellar gravitational collapse trigger PHRASE u. The 
events selected in this (3 range are rejected since their 
pulse width is smaller than the expected counter cross- 
ing time or since the light produced is lower than that 
expected for a monopoles. The flux upper limit at 90 % 
C.L. is 4.3 x 10~ 16 cm- 2 s~ 4 sr^Jpurve "D" in fig. |). 
The technique is fully discussed inllJ. 

Finally, in the high velocity range ((3 > 0.1) the data 
collected by the muon trigger ERPE3 are used. All events 
are rejected since the measured energy deposit in_two 
counter layers is much lower than the energy losscl ex- 
pected for a fast monopole. The 90 % C.L. flux up 
limit is 4.4 x 10 -15 cm~ 2 s~ 4 sr~ 4 (curve "C" in fig. 

3.2 Search with the streamer tubes subdetector 
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The MACRO streamer tubes &EI are filled with a mix- 
ture of He (73 %) and n-pentane (27 %); the Helium 
was chosen to allow the detection of slqw (/3 ;$ 10~ 3 ) 
monopoles by the Drell-Penning effect 113. The hits on 
the streamer tubes system and the charge collected on 
each tube provide measurements oLtrack, velocity and 
energy loss of an ionizing particle til. The data were 
collected from February 1992 to October 1997. The 
monopole analysis is based on the search for clean single 
tracks of well reconstructed velocity; it was checked that 
the trigger selection and analysis procedure are veloc- 
ity independent. No candidated survived; the flux upper 
limit (90 % C.L.) is $ < 4.5 x 10~ 16 cm- 2 s" 1 sr 
1.1 xlO" 4 < < 5xl0- 3 (curve "Streamer" infig.|)t 

3.3 Search with the nuclear track subdetector 

The MACRO nuclear track subdetector is made by three 
layers of LEXAN and three layers of CR39; its total sur- 
face is 1263 m 2 and its acceptance for fast monopoles 
is w 7100 m 2 sr. A calibration of the CR39 with slow 
and fast ions showed that its response depends on the 
restricted energy loss only Q. The track-etch subdetec- 
tor is used in a stand-alone or in a "triggered" mode by 
the streamer tubes and the scintillator systems. A total 
surface of 181 m 2 was etched, with an average exposure 
time of 7.23 years; the flux upper limits (90 % C.L.) arc 
- 0.88 and - 1.3 x 10" 15 cm" 2 s" 1 sr" 1 at /3 ~_L and 
(3 ~ 10~ 4 respectively (curves "CR39" in fig. 

3-4 Combined search 

The fast monopoles are expected to release a huge 
amount of energy by ionization/excitation without pro- 



ducing showers for (3 < 0.99. A monopole search based 
on the energy deposition only can be seriously affected 
by a large background due to showering cosmic rays, but 
such a background is efficiently rejected by a combined 
analysis. This analysis uses at the same time the data 
collected by the scintillators and by the streamer tubes, 
requiring large photomultiplier pulses, an isolated single 
track and a high streamer charge per unit path length. 
Only few (~ 5/year) events survive and are looked for in 
the appropriate track-etch sheets. In 667 days of live time 
no candidates were found; the 90 % C.L. flux upper limit 
is 1.5 x 10~ 15 cm' 2 ^- 1 sr" 1 for 5 x 10~ 3 < f3 < 0.99 
(curve "E" in fig. 

3.5 Conclusions about magnetic monopoles 

We show in fig. [I] the flux upper limits obtained using 
the various MACRO subdetectors. Since each subde- 
tector can rule out, within its acceptance and sensitiv- 
ity, a potential candidate from the others, we obtain a 
global MACRO limit (curve "GLOBAL" in fig. Q), as 
an "OR" combination of the separate results. The pre- 
scriptions used for this combination are described in 13. 
In fig. H we compare the MACRO combined result with 
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Figure 1: Magnetic monopole flux upper limits obtained using the 
various MACRO subdetectors. 



that obtained by other experimentsElli3'Bli!lol23'Eil and 
with the Parker Bound. Our limit is at the level of 0.3 
times the Parker Bound for (3 > 10~ 4 and is the best 
existing for 10~ 4 < (3 < 5 x 10~ 2 . 

4 Nuclearites 

The results obtained using the liquid scintillator and the 
nuclear track subdetectors can be, at least in part, ex- 
trapolated to the search for nuclearites E3, hypothesized 
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Figure 2: Magnetic monopolc flux upper limits obtained by 
MACRO and by other experiments. 



Figure 3: Nuclearite flux upper limits obtained by MACRO and 
by other experiments; the dark matter bound is also shown. 



nuggets of strange quark matter (the streamer tubes are 
not sensitive to_nuclearites because of the low density of 
the filling gas E3) . ._By studying the mechanism of nu- 
clearite energy losscJ it was..shown that the scintillators 
are sensitive to nuclearitesc3 down to /3 ^ 10~ 4 and the 
CR39EI down to fi w 1CT 5 . The MACRO flux upper 
limits are ~4x 10~ 16 cm~ 2 s _1 sr _1 for Mn > 0.1 g 
and about 2 times higher for Mjv < 0.1 g (the lighter 
nuclearites cannot traverse the Earth and then only the 
down-going flux must be considered). Assuming a nu- 
clearite velocity at the ground level j3 — 2 x 10~ 3 wc 
comp ared nu r limit with the results of other experi- 
mentsEj'tSo and with the dark matter bound; our results 
look competitive in a large range of nuclearite masses 
(10 14 GeV/c 2 < Mjy £ 10 21 GeJ^c 2 ) (see fig. ||.) For 
further details on this search seec3. 

5 WIMPs 

The Weakly Interacting Massive Particles (WIMPs) are 
important candidates for the Cold non-baryonic part of 
the Dark Matter in the Universetj. Between the various 
Cold Dark Matter candidates subject to weak interac- 
tions one of -the most promising is the supersymmetric 
neutralino xEj. 

In supersymmetric theories where the R parity is 
conserved there exists a lightest stable supersymmetric 

6 The limit is prolonged above the Dark matter bound 
to show the transition to an isotropic flux for M > 
0.1 g (« 4.5 x 10 22 GeV/c 2 ) 



particle (LSP), which is the natural candidate for the 
Dark Matter since its expected density is close to the 
critical one: FIlsp ~ 1- In many theories the LSP is 
the neutralino x, the simplest linear combination of the 
gaugino and higgsino eigenstates. The \ mass depends 
on the supersymmetric parameters, as, for instance, the 
gaugino and higgsino mass parameters M\, Mi and fi 
and the ratio of the Higgs doublet vacuum expectation 
values tan (3. These parameters are constrained by accel- 
erator searches and a lowerJimit on was set by LEP 
2 data: m% > 20 30 GeVEI The search for WIMPs in 
underground detectors can probe complementary regions 
of the parameter space. 

5.1 WIMP searches m MACRO 

The WIMPs are indirectly searched in MACRO by using 
upward going muons. A WIMP intercepting a celestial 
body can lose its energy and be trapped in the core of 
this body and annihilate with an other WIMP; the de- 
cay of the annihilation products produces high-energy 
i/s which can be detected in an underground detector 
as upward going muons. Then, a statistically significant 
excess of upward going muons from the direction of a 
celestial body can be a hint for a WIMP-WIMP annihi- 
lation in the core of that body. Some possible traps were 
proposed and some calculations of the WIMP annihila- 
tion rate in the Earth and Sun and of the c^pcsppnding 
upward going muon fluxes were performed oEjlj. This 
indirect search achieves a better signal to noise ratio for 
high WIMP masses, since the higher the WIMP mass, 
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the more the upgoing muon follows the parent neutrino 
direction. j-IheJtechnique was already used by other ex- 
periments&OEJ. 

The upward going muons in MACRO are selected 
by the time-of-flight technique, requiring a track in the 
streamer tubes and a time of flight between the scintilla- 
tion counter layers consistent with a j3 ~ 1 particle from 
below. Details on the selection criteria can be found inE3. 

5.2 Search for WIMP annihilation in the Earth 

Upward going muons from WIMP annihilation in the 
Earth core are searched in angular cones (3 -j- 30° wide) 
around the vertical. We used 517 events, collected in 
3.1 years of live time, requiring a minimum crossing of 
200 g cm~ 2 of rock absorber. After background subtrac- 
tion the number of selected events is 487 ± 22 sta t, to be 
compared with 653 ± lllt/ leor expected from a Monte 
Carlo calculation. Most of the deficit lies around the ver- 
tical, where the WIMP annihilation signal is expected, 
but also where the efficiency and acceptance of the appa- 
ratus are best known. To set a conservative limit on the 
WIMP flux we assumed that the number of measured and 
expected events are equal and normalized the expected 
distribution to the factor 0.85, corresponding to the ra- 
tio between measured and expected events for 9 > 30°. 
With this prescription our limits on the WIMP flux from 
the Earth range from 0.4 to 2.3 x 10~ 14 cm~ 2 s _1 for 
angular windows from 3 to 30°. 



5.3 Search for WIMP annihilation in the Sun 

In the search for WIMPs from the Sun, since the back- 
ground for moving sources is lower than for steady 
sources, we used an enlarged sample (762 events) which 
includes also muons partially contained in the apparatus 
(produced by neutrino interactions in the absorber in the 
MACRO lower part) and muons crossing < 200 g cm -2 of 
rock absorber. The simulation was obtained by the data 
themselves to include properly the effects of the partially 
contained events and the arrival times were extracted 
randomly during the whole measurement time to take 
into account possible drifts of detection efficiency. The 
angular distribution of the upward going muon events do 
not show any significant excess around the Sun direction. 
Then, we set upper limits on the WIMP flux from the 
Sun ranging from 1.7 to 6.0 x 10~ 14 cm~ 2 s _1 in the 
angular window 3 -j- 30°. 



5.4 Limits on the WIMP fluxes and \ mass from the 
angular distributions 

As already stated, the angle between the upward going 
muon and the neutrino from the WIMP annihilation de- 
pends mainly on the WIMP mass. We performed a full 



Monte Carlo calculation of the expected angle between 
the upward going muon and the Earth or Sun directions 
for neutralino masses from 60 GeV to 1000 GeV. This 
calculation uses the shape of the upward going muon sig- 
nals from x pnX annihilation in the Earth and the Sun 
computed in Ed, propagates the muons through thfi rock 
to the apparatus (with the cross sections given inEZI) and 
includes the angular smearing produced by the experi- 
mental resolution. As expected, the angle becomes nar- 
rower when increases, but also in the less favourable 
case more than 90 % of the signal is contained in an an- 
gular cone 9 < 15° around the selected direction. For 
each m-^ we determine the cone which collects the 90 % 
of the expected signal and define a corresponding 90 % 
C.L. limit on the upward going muon flux from X — X 
annihilation. These limits for the Sun are shcpsui in fig. ^ 
and compared with the fluxes computed in Ell, varying 
some supersymmetric parameters (Mi, fj,, tan/3 etc.). 
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Figure 4: Upward going muon flux from the Sun vs m-j computed 
in 31 . The solid line is the 90% C.L. MACRO flux limit. 

Fig. ^ shows that, thanks to the improved statistics and 
exposure, our data.s±art to constrain the supersymmetric 
theoretical models cB. 

6 LIPs 

Fractionally charged particles have been actively 
searched for since many years, but without success; the 
detection of such particles would be a proof of their ex- 
istence and/or of a lack of the confinement hypothesis 
under some circumstances. The family of the possible 
fractionally charged particles includes the quarksEJ (with 
charge |e| /3 and 2/3 |e|) and a variety of other particles 
predicted by Grand Unified Theories, with charges rang- 
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ing from 1/5 |e|upto2/3 |e|00@. Particles with frac- 
tional charge deposit less energy than particle with unit 
charge since the energy loss is proportional to the square 
of the charge; such particles are called Lightly Ionizing 
Particles (LIPs). 

A specialized trigger was developed in MACRO for 
the LIP search; this trigger uses the low energy events 
collected by the PHRASE system and performs a four- 
fold coincidence between the signals coming from three 
counters (each one in a different scintillator layer) and 
from the streamer tubes. The LIP trigger is sensitive 
down to \Q/e\ = 1/5; the corresponding energy loss in the 
MACRO counters is AE « 1.6 MeV. The LIP trigger 
gives the stop to the WFD system; the WFD data pro- 
vide high quality measurements of the energy and timing 
of the events. For each LIP event, the maximum en- 
ergy loss in the three scintillator layers is computed and 
compared with that expected for a fractionally charged 
particle. The efficiency as a function of the charge is 
obtained by a Monte Carlo simulation which takes into 
account the trigger behaviour close to the threshold and 
the cosmic muon background. No candidate survived out 
of 1.2 million triggers; the 90% C. L. upper limit for an 
isotropic flux is <I> < 9.2 x 1CT 15 cm _2 s _1 sr _1 . This 
limit is shown in fig. (solid line-Land compared with 
that set by other collaborationsBo. 
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Figure 5: 90 % C.L. flux upper limits for LIP set by MACRO (solid 
line) and by other experiments as a function of the LIP charge. 

We stress the fact that our experiment is the first 
one to be sensitive down to \Q/e\ = 1/5 and that we 
can perform a high quality search thanks to the com- 
bined signature provided by the streamer tubes and the 
scintillators: only few events over 1.2 millions needed a 
hand scanning (instead of one event over few thousands 
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